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Abstract — In this paper it is shown how one can use Bessel beams to obtain a 
stationary localized wavefield with high transverse localization, and whose longitudinal 
intensity pattern can assume any desired shape within a chosen interval < z < L of 
the propagation axis. This intensity envelope remains static, i.e., with velocity v — 0; 
and because of this we call "Frozen Waves" such news solutions to the wave equations 
(and, in particular, to the Maxwell equations). These solutions can be used in many 
different and interesting applications, as optical tweezers, atom guides, optical or acoustic 
bistouries, various important medical purposes, etc. 
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1. — Introduction 

Since many years a theory of localized waves (LW), or nondiffracting waves, has been 
developed, generalized, and experimentally verified in many fields as optics, microwaves 
and acoustics[l]. These waves have the surprising characteristics of resisting the diffraction 
effects for long distances, i.e., of possessing a large depth of field. 

These waves can be divided into two classes, the localized beams, and the localized 

pulses. With regard to the beams, the most popular is the Bessel beam[l]. 

(t) Work supported by FAPESP (Brazil), ; previously available as e-print E-mail address 

for contacts: mzamboni@dmo.fee.unicamp.br 
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Much work have been made about the properties and applications of single Bessel 
beams. By contrast, only a few papers have been addressed to the properties and ap- 
plications of superpositions of Bessel beams with the same frequency, but with different 
longitudunal wave numbers. The few works on this subject have shown some surprising 
possibilities related with this type of superpositions, mainly the possibility of controlling 
the transverse shape of the resulting beam [2, 3]. The other important point, i.e., that of 
controlling the longitudinal shape, has been very rarely analyzed, and the relevant papers 
have been confined to numerical optimization processes[4,5] to find out one appropriate 
computer-generated hologram. 

In this work we develop a very simple method** that makes possible the control of 
the beam intensity longitudinal shape within a chosen interval < z < L, where z is the 
propagation axis and L can be much greater than the wavelength A of the monochromatic 
light which is being used. Inside such a space interval, we can construct a stationary 
envelope with many different shapes, including one or more high-intensity peaks (with 
distances between them much larger than A). This intensity envelope remains static, i.e., 
with velocity V = 0; and because of this we call "Frozen Waves" such news solutions to 
the wave equations (and, in particular, to the Maxwell equations). 

We also suggest a simple apparatus capable of generating these stationary fields. 

Static wave solutions like these can have many different and interesting applications, 
as optical tweezers, atom guides, optical or acoustic bisturies, electromagnetic or ultra- 
sound high-intensity fields for various important medical purposes, etc..** 

2. — The mathematical methodology ** 

We start with the well known axis-symmetric Bessel beam 



i>(p,z,t) = J (k p p)e i(3z e 



(1) 



with 



uj 2 

k p = — - , 



(2) 



** Patent pending. 
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where u, k p and (3 are the angular frequency, the transverse and the longitudinal wave 
numbers, respectively. We also impose the conditions 



u//3 > and k 2 p > (3) 

to ensure forward propagation only, as well as a physical behavior of the Bessel function. 

Now, let us make a superposition of 2N + 1 Bessel beams with the same frequency 
u> , but with different (and still unknown) longitudinal wave numbers (3 n : 

N 

*(p,z,t) = e~^ ot A n Uk pn p)e l ^\ (4) 

n=—N 

where A n are constant coefficients. For each n, the parameters uq, k pn and (5 n must satisfy 
Eq.(J2J), and, because of conditions (jSJ), when considering u > 0, we must have 

0<Pn<~ (5) 

c 

Now our goal is to find out the values of the longitudinal wave numbers f3 n and of the 
coefficients A n in order to reproduce approximately, inside the interval < z < L (on the 
axis p = 0), a chosen longitudinal intensity pattern that we call |F(z)| 2 . In other words, 
we want to have 

N 

A n e iPnZ ^F(z) with 0<z<L (6) 

n=—N 

Following Eq.©, one might be tempted to take (3 n = 2-nn/L, thus obtaining a trun- 
cated Fourier series, expected to represent the desired pattern F(z). Superpositions of 
Bessel beams with (3 n = 2im/L has been actually used in some works to obtain a large 
set of transverse amplitude profiles[2]. However, for our purposes, this choice is not ap- 
propriate due two principal reasons: 1) It yields negative values for f3 n (when n < 0), 
which implies backwards propagating components (since ujq > 0); 2) In the cases when 
L >> Ao, which are of our interest here, the main terms of the series would correspond 
to very small values of (3 n , which results in a very short field depth of the corresponding 
Bessel beams(when generated by finite apertures), impeding the creation of the desired 
envelopes far form the source. 
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Therefore, we need to make a better choice for the values of /3 n , which allows forward 
propagation components only, and a good depth of field. This problem can be solved by 
putting 

(3 n = Q + ^n, (7) 

where Q > is a value to be chosen (as we shall see) according to the given experimental 
situation, and the desired degree of transverse field localization. Due to Eq.flHJ), we get 

0<Q±^N<^ (8) 

1j c 

Inequation (jSJ) determines the maximum value of n, that we call N, once Q, L and t^o 
have been chosen. 

As a consequence, for getting a longitudinal intensity pattern approximately equal to 
the desired one, F(z), in the interval < z < L, Eq.(@J) should be rewritten as: 

N 

V(p = 0,z,t) = e - iuJot e iQz J2 A n e i2 -? nz , (9) 

n=-N 

with 

A n = 1 [ L F(z)e- i ^ nz dz (10) 
LJo 

Obviously, one obtains only an approximation to the desired longitudinal pattern, 
because the trigonometric series Q is necessarily truncated. Its total number of terms, 
let us repeat, will be fixed once the values of Q, L and u are chosen. 

When p ^ 0, the wave field ^(p, z, t) becomes 

*(p,z,t) = e -^ l e iQz J2 A n J {k pnP )e^ nz , (11) 

n=-N 

with 

Kn = -[Q+ —J (12) 

The coefficients A n will yield the amplitudes and the relative phases of each Bessel 
beam in the superposition. 
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Because we are adding together zero order Bessel functions, we can expect a high field 
concentration around p = 0. 



3. — Some examples 

In this section we shall present two examples of our methodology. 

Let us suppose that we want an optical wave field with Ao = 0.632 /im, that is, with 
Uq = 2.98 10 15 Hz), whose longitudinal pattern (along its 2-axis) in the range < z < L 
is given by the function 



{ (z-jiKz-ja) 



F(z) 



-4 



(z - l 5 )(z - k) 
(k - h) 2 





for l\ < z < l 2 
in l 3 < z < I4 
for Z 5 < z < Iq 
elsewhere , 



(13) 



where h = L/10, l 2 = 3L/10, h = 4L/10, l 4 = 6L/10, h = 7L/10 and l 6 = 9L/10. 
In other words, the desired longitudinal shape, in the range < z < L, is a parabolic 
function for li < z < I2, a unitary step funtion for Z3 < z < U, and again a parabola 
in the interval h < z < Iq, it being zero elsewhere (in the interval < z < L). In this 
example, let us put L = 0.5 m. 

We can then calculate the coefficients A n , which appear in the superposition (JUJ), by 
inserting Eq. (fT3"j) into Eq.ffTUj). Let us choose, for instance, Q = 0.9998 u /c: This choice 
allows the maximum value iV = 158 of n, as one can infer from Eq.(jSJ). Let us specify 
that, in such a case, one is not obliged to use just N = 158, but one can adopt for iV any 
values smaller than it; more in general, any value smaller than that calculated via Eq. ({%)). 
Of course, on using the maximum value allowed for N, one will get a better result. 

In the present case, let us adopt the value iV = 20. In Fig. 1(a) we compare the 
intensity of the desired longitudinal function F(z) with that of the Frozen Wave (FW), 
^(p = 0, z, t), obtained from Eq. (JUJ) by using the mentioned value N = 20. 

One can verify that a good agreement between the desired longitudinal behaviour and 
our approximate Frozen Wave is already obtained with iV = 20. Obviously, the use of 
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(a) 



(b) 



Figure 1: (a) Comparison between the intensity of the desired longitudinal function F(z) 
and that of our Frozen Wave (FW), \I/(p = 0,z,t), obtained from Eq.Q. The solid 
line represents the function F(z), and the dotted one our FW. (b) 3D-plot of the field 
intensity of the FW chosen in this case by us. 



higher values for N will improve the approximation. 

Fig. 1(b) shows the 3D- intensity of our FW, given by Eq. (jllj) . One can observe that 
this field possesses the desired longitudinal pattern, while being endowed with a good 
transverse localization. 

We can expect that, for a desired longitudinal pattern of the field intensity, on choosing 
smaller values of the parameter Q one will get FWs with higher transverse width (for the 
same number of terms in the series ([lip), because of the fact that the Bessel beams in 
(lllj) will possess a larger transverse wave number, and consequently higher transverse 
concentrations. We can verify this expectation on considering, for instance, a desired 
longitudinal pattern, in the range < z < L, given by the function 



(z-h)(z-l 2 ) 

F(z) = { ^~ ^ , (14) 

in the otherwise 

with 11 = L/2 — AL and 12 = L/2 + AL. Such a function has a parabolic shape, with 



6 



the peak centered at L/2 and a width of 2AL. By adopting A = 0.632 /zm (that is, 
ojq = 2.98 10 15 Hz), let us use the superposition (JTTJ) with two different values of Q: we 
shall obtain two different FWs that, in spite of having the same longitudinal intensity 
pattern, will have different transverse localizations. Namely, let us consider L = 0.5 m 
and AL = L/50, and the two values Q = 0.99996u;o/c and Q = 0.99980cj /c. In both 
cases the coefficients A n will be the same, calculated from Eq. fjlOJ) . on using this time the 
value N = 30 in the superposition (jf fj) . The results are shown in Figures (2a) and (2b). 
One can observe that both FWs have the (same) longitudinal intensity pattern, but the 
one with the smaller Q is endowed with the higher transverse localization. 



m 2 




Figure 2: (a) The Frozen Wave with Q = 0.99996co>o/c and N = 30, approximately 
reproducing the chosen longitudinal pattern represented by Eq.(^3J). (b) A different 
Frozen wave, now with Q = 0.99980u;o/c (but still with iV = 30) forwarding the same 
longitudinal pattern. We can observe that in this case (with a lower value for Q) a higher 
transverse localization is obtained. 



4. — Generation of Frozen Waves 

Concerning the generation of Frozen Waves, we have to recall that the superpositions 
which define them, consists of sums of Bessel beams. Let us also recall that a Bessel 
beam, when generated by finite apertures (as it must be, in any real situations), maintains 
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its nondiffr acting properties till a certain distance only (its field depth), given by 



^ = -^tt , (15) 
tanfl v ; 

where R is the aperture radius and 9 is the so-called axicon angle, related with the 

longitudinal wave number by the known expressionfl] cos# = c(3/u. 

So, given an apparatus whatsoever capable of generating a single (truncated) Bessel 

beam, we can use an array of such apparatuses to generate a sum of them, with the 

appropriate longitudinal wave numbers and amplitudes/phases (as required by Eq.([TT|)). 

thus producing the desired FW. Here, it is worthwhile to notice that we shall be able 

to generate the desired FW in the the range < z < L if all Bessel beams entering the 

superposition (JTTJl are able to reach the distance L resisting the diffraction effects. We 

can guarantee this if L < Z min , where Z min is the field depth of the Bessel beam with the 

smallest longitudinal wave number /3 n =-N — Q — 2nN/L, that is, with the shortest depth 

of field. In such a way, once we have the values of L, ujq, Q, N, from Eq. (jl5|) and the 

above considerations it results that the radius R of the finite aperture has to be 



R> L. 



""I (16) 



^n=-N 



\ c 2 A 

The simplest apparatus capable of generating a Bessel beam is that adopted by Durnin 
et al. [6], which consists in an annular slit located at the focus of a convergent lens and 
illuminated by a cw laser. Then, an array of such annular rings, with the appropriate 
radii and transfer functions able to yield both the correct longitudinal wave numbers* and 
the coefficients A n of the fundamental superposition (jllj) . can generate the desired FW. 
This questions will be analyze in more detail elsewhere. 

Obviously, other powerful tools, like the computer generated holograms (ROACH's 
approach, for instance), may be used to generated our FWs. 

5. — Conclusions 

In this work we have shown how Bessel beams can be used to obtain stationary lo- 
calized wave fields, with high transverse localization, whose longitudinal intensity pattern 



*Once a value for Q has been chosen. 
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can assume any desired shape within a chosen space interval < z < L. The produced 
envelope remains static, i.e., with velocity V = 0, and because of this we have called 
Frozen Waves such news solutions. 

The present results can find applications in many fields:** For instance, in the optical 
tweezers modelling, since we can construct stationary optical fields with a great variety 
of shapes, capable, e.g., of trapping particles or tiny objects at different locations. This 
topic is being studied and will be reported elsewhere. 
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